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Abstract 



Given a random walk (5„)„gz defined for a doubly infinite sequence of times, we let the time 
parameter {nk)keN itself be a process with values in Z and call {Snk)ken random walk at ran- 
dom time. We show that under suitable conditions, it scales to an (ff-sssi)-time Levy motion, a 
generaUzation of iterated Brownian motion. 

In Khoshnevisan and Lewis (1999), a normalized variation result for iterated Brownian motion 
was stated which alternated between even and odd orders of variation. This result "suggested the 
existence of a form of measure-theoretic duality" between iterated Brownian motion and a Brownian 
motion in random scenery. We show that a random walk at random time can be considered a random 
walk in "alternating" scenery, thus hinting at a mechanism behind this duahty. 

Following Cohen and Samorodnitsky (2006), we also consider alternating random reward schema 
(random reward schema are certain sums of random walks in random scenery) associated to random 
walks at random times. Random reward schema scale to local time fractional stable motions which 
are ff-sssi symmetric a-stable processes with self-similarity exponents H > 1/a. We show that 
the altemating random reward schema scale to indicator fractional stable motions which are ff-sssi 
symmetric a-stable processes with self- similarity exponents H < 1/a. 

We also show that one may use a recursion on the subordinating random time processes to get 
new local time and indicator fractional stable motions and new stable processes in random scenery 
or at random times. When a = 2, the fractional stable motions given by the recursion are fractional 
Brownian motions with dyadic H G (0, 1). Finally, we show that "undoing a subordination" via 
a time -change allows one to, in some sense, extract Brownian motion from fractional Brownian 
motions with H < 1/2. 

Keywords: fractional Brownian motion; random walk in random scenery; random reward schema; 
local time fractional stable motion; self-similar process; iterated process 
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1 Introduction 

Let (<), i = 1, 2, 3 be three independent Brownian motions, and let a two-sided Brownian motion 
be defined by 



motion (IBM). It can be thought of as a two-sided Brownian motion which is non-monotonically 
"subordinated" to another Brownian motion. This process was also used by IIDM921 to study the 
Bahadur- Kiefer process. Also, a variant of IBM, where the pure imaginary process ii?'^^^ (~t) was 
substituted for t < 0, was utilized by l|Fun79l to study the PDE: 

du 1 d^u 

Recently, more general processes at random times called a-time Brownian motions and a-time 
fractional Brownian motions were introduced in IINan06l INWX 111. In these works (along with sev- 
eral references therein), the connection between processes at random times and various PDE's was 
studied, along with the local time and path properties of the iterated processes. In a different direc- 
tion, the scaling and asymptotic density of a discretized version of IBM called iterated random walk 
was analyzed in the physics literature IITur04L 

In this work, we consider generalizations of the iterated random walk which we call random 
walks at random times (RWRT) and dependent walks at random times (DWRT) and relate them with 
a different portion of the probability literature concerning random walks in random scenery. This 
relation was first noted in lKL96i who stated that there was "a surprising connection between the 
variations [of IBM] and H. Kesten and F. Spitzer's Brownian motion in random scenery." Later, in 
IIKL99I , a form of measure-theoretic duality was shown between the two processes. Here, we present 
a mechanism on the discrete level which shows a connection between the two processes. 

We show that under suitable conditions, the scaling limits of RWRT and DWRT are (_ff-sssi)- 
time a-stable Levy motions, a new class of processes at random times. If X{t) is a two-sided 
a-stable Levy motion defined similarly to ([T]i, and Yt is an independent a-stable Levy motion then 
we call X{Yt) an iterated Levy motion. If, more generally, Yt is an independent i7-self-similar, 
stationary-increment process (sssi), then an (77-sssi)-time a-stable Levy motion is given by X{Yt). 
Assuming < 7? < 1, we will see that X{Yt) is an i7/a-sssi process with Hurst exponent less 
than 1/a. They naturally complement stable processes in random scenery which are the limiting 
continuous processes of IIKS79I and MWan03l and which have Hurst exponents greater than 1/a 
( IIWan 03 1 considered only the case a = 2, but this was extended to a < 2 by iCD09l ). 

Random walks in random scenery (RWRS) and their scaling limits, stable processes in random 
scenery, were first introduced independently in BKS79I |Bor79l . The purpose of IIKS79I was to in- 
troduce a new class of sssi processes given by the scaling limits of RWRS. The scaling limits have 
integral representations as stable integrals of local time kernels (of a process Yt). When the ran- 
dom scenery are a-stable laws, they scale to the a-stable random measure against which the local 
time kernel is integrated. In comparison, there is also an integral representation of (i/-sssi)-time a- 
stable Levy motions given by the stable integration of random kernels of type l[o, Yt] against a-stable 
random measures. 

When Yt is a generic 7?-sssi process, the stable processes in random scenery discussed above 
also include the model of IIWan03l . Wang used "dependent walks" to collect the scenery, instead of 
random walks, leading to a dependent walk in random scenery (DWRS). In particular, the dependent 
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In BBur92| . Burdzy studied the process 
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Figure 1: a S (1, 2) : LT = local time, I = indicator, L = linear, RH = real harmonizable 



walks he used were discrete-time Gaussian processes known to scale to fractional Brownian motion 
(fBm). 

Random reward schema are sums of independent copies of discrete processes in random scenery. 
In aCS06l iDGPOSl ICD09I it was shown that the random reward schema of RWRS and DWRS scale 
to iJ-sssi symmetric a-stable (SaS) processes called local time fractional SaS motions (with H > 
1/a). In this work, we show that the scaling limits of random reward schema for RWRT and DWRT 
are ff-sssi SaS processes called indicator fractional SaS motions (with H < 1/a) which were 
introduced in I Junl 111 . 

Note that fBm is the only sssi Gaussian process. Thus, when the scenery has finite variance and 
a — 2, local time fractional SaS motions and indicator fractional SaS motions reduce to fBm with 
H > 1/2 and H < 1/2, respectively. 

As will be seen in Section |2] the mechanism behind the connection between local time and 
indicator fractional stable motions is the same as the mechanism which connects Brownian motion in 
random scenery (BMRS) with IBM. In effect, the mechanism shows that the indicator kernels of the 
latter processes can be thought of as "alternating" versions of the local time kernels of the former 

Together, local time fractional SaS motions and indicator fractional stable motions form a class 
of fractional stable motions (i7-sssi SaS processes) which may be thought of as one of several 
generalizations of fractional Brownian motion. Their increment processes are stationary and have the 
ergodic -theoretic property of being null conservative, a concept introduced in IISam05l . This property 
distinguishes them from fractional stable motions which have dissipative or positive conservative 
increment processes. The most well-known examples of fractional stable motions with dissipative 
or positive conservative increment processes are the linear fractional stable motions and the real 
harmonizable stable motions, respectively, as can be seen in the figure below. 

We also consider single-scenery random reward schema introduced in IIDGP09L Here we again 
take sums of identically distributed RWRTs or DWRTs. However the copies have a dependence 
structure since they use the same "single scenery". This dependence will be made more explicit 
below. The scaling limits of single-scenery random reward schema of RWRS and DWRS no longer 
have stationary increments, however, they are easily seen to be H-ss SaS processes with H > 1/a. 
Similarly, the scaling limits of single-scenery random reward schema of RWRT and DWRT are H-ss 
SaS processes with H < 1/a. 

Finally, we also present a recursive construction of some local time and indicator fractional stable 
motions. In particular, we show that at each step of the recursion, the local times exist and are in 

{^l X M.). The recursively defined processes give the first examples of local time fractional stable 
motions for which the processes collecting the scenery are neither fBm nor /3-stable Levy motions. 
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In the case a = 2, the processes are given by integrals against Gaussian random measures, and the 
recursion constructs fBm, of any dyadic Hurst parameter, using one Brownian motion and a countable 
family of independent random Gaussian measures. 

As mentioned above, RWRT and, in particular, its scaling limit are in some sense non-monotonically 
subordinated processes. Usually one may not undo a subordination- for example one can embed a 
stable process in Brownian motion, but cannot extract Brownian motion from the stable process since 
the filtration is strictly smaller. However we will see that when the scaling limit of the random time 
process, Yt, is fBm, one can undo the subordination using the time-change Ts — inft>o{< : Yt — s}. 
Extending such a time-change procedure to the kernels of indicator fractional stable motions when 
a = 2, we find that one can, in some sense, extract Brownian motion from fractional Brownian 
motions satisfying H < 1/2. 

The rest of the paper is arranged as follows. In Section|2]we describe RWRTs and RWRSs. We 
also describe their respective random reward schema and scaling limits. The section ends with a 
statement describing new scaling limit results. The proofs of the scaling weak convergence results 
are given in Section [3] In Section |4] we describe the recursive construction mentioned above, and 
complete the nontrivial task of showing that the recursion produces processes that are well-defined. 
The main component of this task is showing that the local times exist and are in (£7 x M) . Finally, in 
Section|5]we explain how to extract Brownian motion from fBm with any Hurst parameter satisfying 
H < 1/2. 

2 Discrete and continuous models 

2.1 Random walks at random times and alternating random reward schema 

We start with a simple description of RWRS. Let {i]a{k)}kGZ be a set of i.i.d. symmetric random 
variables in the domain of attraction of an SaS law, a G (0, 2] with scale parameter a = 1. The 
family {r]a{k)} depicts the scenery associated to the vertices of Z. Let 

n 

W{n):=Y,^0{k) (3) 

fe=i 

be a symmetric random walk on Z with steps £,i3{k) in the domain of attraction of an S/3S law, 
P e (1. 2]. The random walk roams amidst the scenery {i]a{k)} which are independent from the 
steps {Cpik)}. 

The cumulative scenery process 

n 

Zn = Z.n{Va,W) ^ (M^(fc)) (4) 

is called a random walk in random scenery. The scenery {ria{k)} can alternatively be thought of as 
random rewards collected by the random walk when it visits vertex k. 

We note that some authors call the pair {W,ria{W)) a RWRS process (for example lldHS06ll ). 
Since most of the papers cited in this work refer to Q as the RWRS, we stick with this notation. 

Wang IIWan03l considered a slight modification of RWRS by using a discrete approximation of 
a Gaussian process instead of a random walk: 



n 
k=l 



(5) 
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Here [•] is the ceiling function and G/f (fc) is the partial sum of a stationary Gaussian process Xk 
with correlations r{j — k) = EXjXk satisfying 

n n 

^^r(j-fc)^n2^, (6) 

3 = 1 k=l 

where < < 1. In addition to (|5]l, there have been myriad generalizations of Q and we refer the 
reader to the introduction of IGPPIOII for a nice summary of such generalizations. 

We refer to (|5]) as a dependent walk in random scenery (DWRS). In general, we consider Z„ {rja , Wh) 
for which the collecting process Wr in) has stationary increments and also satisfies the following 
scaling limit properties: 

(z) \imn^^n-"WH{\nt\)^Yt in 2?([0, c»)) , 
(SLP) I (ii) Yt is a non-degenerate iJ-sssi process (Fq = by self-similarity), 
(m) ElYfl < oo, 

where I?([0, oo)) is equipped with the usual Skorohod topology (also called the Ji -topology). 

The condition that Yt be sssi guarantees that Z„ scales to an sssi process as well, and this was 
in fact the original motivation of introducing Z„ in [KS791 . Note that we use the stable parameter 
a S (0, 2] for the scenery/rewards and consequently the increments of the RWRS/DWRS, however 
we reserve the stable parameter /3 S (1, 2] for the increments of the collecting process (note that we 
require /3 > 1 in order to guarantee E|Yf | < oo). 

We introduce a variant of Z„ in which the rewards alternate in sign and are associated with 
edges instead of vertices. In our variant of RWRS, we use symmetric rewards {?/Q(e)} together with 
signs {(Te}, G {—1, +1} associatcd to the edge set of Z. At time zero, all signs are plus one, 
(Te(0) — +1, however ((Te(n))„>o is a process determined by the collecting process in a manner 
discussed below. 

Consider a discrete collecting process Wnin) satisfying condition (SLP). Note that our definition 
allows |tV//(rt) — WH{n — 1)| to be greater than one. Let be the set of connected edges traversed 
on the nth step of Wnin), i.e. the set of edges between Wnin — 1) and WHin) (thus f„ has 
cardinality |Wi/(ri) — Wnin — 1)|). At the nth step, the process Wnin) 

• earns the signed rewards CTe(n — 1) • ria{e) of all edges e £ and then 

• reverses the sign Ce of each e € so that it will receive the exact opposite reward the next 
time it traverses e. 

A (dependent) random walk at random time (DWRT/RWRT) with a non-monotonically subordi- 
nating random time process Wh (n) is a process 

n 

An = An{Va,WH) ■= ^ '^•^(^ " ^) ' '?"(^)- 

k—1 e^£-k 

where ae{k) E {—1, +1} is the sign of e at time k. 

To explain the name of the process, consider that in an RWRT, due to cancelation, each reward 
rjaie) contributes either one or zero net terms to the sum (|7|. When e is to the right of the origin, the 
number of net terms is one if and only if Wh (n) is to the right of e, and when e is to the left of the 
origin, the number of net terms is one if and only if W^r(n) is to the left of e. It follows that 

An^ ?7a(e), (8) 

ee[0,H'H(n)] 
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where e € [0, x] means that e hes between and x regardless of the sign of x. The partial sum of 
rewards J2ee[o n] Va{s) is just a random walk Sa{n)- If we let Sa{0) — and extend the random 
walk to negative times in the natural way, then thinking of time being determined by the location of 
Wnin), we have 

An = Sc^iWHin)). (9) 

As an aside, if we take (|8]l as our initial definition rather than (j7]i then the rewards may equally well 
be placed on the vertices instead of the edges. The reader may therefore choose to visualize this 
process in any of several ways according to his or her own aesthetic preference. 

The relationship between Z„ and An should be clear. In particular, when the collecting process 
is a simple random walk W{n), a relation is made by using a bijection which assigns to each vertex 
k either the edge lying to its left whenever the previous step of W{n) was in the positive direction 
(right), or the edge lying to its right whenever the previous step of W{n) was in the negative direction 
(left). To extend the relation to other random walks, one must use a modified version of Z„ which, 
when going from a; to y on the nth step, collects a reward not only from y, but all vertices between x 
and y. In view of this relationship between Z„ and An, if Ps is the measure for the random scenery 
and P' is the measure for Wh, then the processes Z„ and An can be defined on the same product 
space with measure Ps x P'. 

There is a further relationship between Z„ and the variations of An which mirrors the connection 
between BMRS and the variations of IBM as presented in IIKL99L In order to explain this relation- 
ship, it will be convenient to let the collecting walk W{n) be a simple random walk and to have the 
rewards for both Z„ and An be attached to the edges of Z, rather than to the vertices. For p g N let 
the p-th variation of An be defined as 



n 



=1 



Theorem 2.1. Suppose the i.i.d. rewards {ijaie)} ore symmetric and have finite pth moments. If p 
is odd, then vlf^ is another RWRT, while ifp is even, then V^^ - nE[<] is a RWRS. In both cases, 
the rewards collected by the processes are given by {^^^^ (e)} where 

Proof. For z > 1 we let 8i denote the edge between W{i — 1) and W{i). We then have 

71 

Zn^^Va{£i), An= ^ Va{e)- (H) 
1=1 ee[0,W{n)] 

Note that 

{A, - = - l)Va{e)f l{£.=e}. (12) 

If p = 2q is even, then the sign {ae{i — l))^"^ in ( 



12 1 is irrelevant. Therefore, 



n n 



Comparing with (111 shows this to be a RWRS with rewards given by {C*^^^ (e)} 

On the other hand, if p = 2(7 + 1 is odd, then the sign {ae{k) — 1))^''^^ — <Je(k) {i — 1) in (|l2| 



causes the same cancelation as we have with RWRT, and since rja is symmetric, there is no longer a 



need to subtract the expectation. Thus, (12 1 yields 
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1=1 ee[O.WH(n)] 



Comparing again with (111 shows this to be a RWRT with rewards given by {C^p) (e)}. □ 



We now compare this with the results of IIKL99L Let Ig denote an IBM, fix an interval [0, t], and 

let 

V^^Ht) = Y.^I{n+l,n) - m^n)r, (15) 

fe=l 

where {Tfe.„ : 1 < fc < 2"i} is an induced random partition of the interval [0, t] (see IIKL991 Sec. 1] 
for details). Among other things, Khoshnevisan and Lewis showed that, when properly renormalized, 
converges in distribution to IBM when p is odd and BMRS when p is even; see Theorems 
3.2, 4.4, 4.5, and the discussion in the middle of page 631. If we consider the natural association 
between BMRS and RWRS on the one hand and between IBM and RWRT on the other, we see that 
the simple Theorem |2. 1 [ provides an intuitive backdrop for the much more difficult results concerning 
the continuous case in |KL99|. 

We now return to study of An in the general case. We will need processes extended to non-integer 
times, and we will therefore denote the linear interpolation of A^ as 

At = A^t\+{t-[mA^t^-A^t\)■ (16) 

Let us now describe the two different random reward schema we will use. Let us start with an 
alternating version of the random reward schema introduced in [ICS06L Let {(W^^(n))„>o}igN be 
independent copies of Wh {n) which are also independent from independent copies of the rewards 
{{?7q ^ (e)}egz}ieN- If (cn) is a sequence of integers such that c„ -> oo, then 

Y^Mv^^^^h) (17) 

i=l 

is an alternating random reward scheme. 

If we instead follow the single-scenery schema of IIDGP09I and use the same single copy of 
rewards {ri'a\e)}e& for each copy of PF^'^(n),then 

YMv^^^^h) (18) 

1=1 

is a single scenery alternating random reward scheme. 

2.2 Scaling limits of random reward schema 

In this section we state some known results concerning the scalings of RWRS and DWRS to sta- 
ble integral representations. These will motivate our results concerning the scalings of RWRT and 
DWRT. 

Let us first recall an important definition. Suppose m is a cr-finite measure on a measurable space 
[E, B), and that 

Bo = e S ; m(A) < c5o}. 
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Definition 2.2. A SaS random measure M with control measure m is a a-additive set function on 
Bq such that for all Ai e So 

1. M(v4i) 

2. M{Ai) and M{A2) are independent whenever Aif] A2 — ^ 

where Sa (c) is an SaS random variable. 
In particular, if f (z L"{E, B, m) then 

[ f{x)M{dx)^S^{\\f{x)U^). (19) 

J E 

Section 3.3 of fST94l contains an introduction to this topic. The immediate importance to us is 
that the scaling limits of RWRS and DWRS are integrals with respect to stable random measures, 
where the integral kernel is the local time of a properly scaled collecting process Wh' (linearly 
interpolated) which is either Gh' or Sp with /3 G (1,2]. The process ^^Wh' {nt) converges weakly 
to a scaling limit, denoted by Yt, which is respectively fBm-iJ' in C([0,oo)) or a /3-stable Levy 
motion in I?([0, 00)). Let (fi', F\ P') be the probability space of Yt- It is known that Yt has a jointly 
continuous local time ^y-(t, x); this was shown for /3-stable Levy motions in | |Boy64[ and for fBm in 
IIBer74i . Moreover, for alH > and all a e (0, 2], Ft satisfies 

E' ( \iy{t,x)\°'dx < 00 (20) 
Jr 

by Theorem 3.1 in |'CS061 and Lemma 2.1 in P DGPOSI . Here we interpret iyW '^^e increasing 
family of random functions which satisfy the occupation time formula 

/ lA{t)ds= f ly{t,X)dX, (21) 
JO J A 

for any Borel set A. 

Let Mo (da;) be an SaS random measure with Lebesgue control measure which is independent 
from Yt- Throughout this subsection we will let 

H = l-H' + H'/a. (22) 

A stable process in random scenery is an _ff-sssi SaS process given by 

Af{Mo,Y) / iyit,x)Mo{dx), t > 0, (23) 

JR 

which is well-defined by ( |20j l (see Chapter 3 of IIST94I ). Recall that ria{k) is in the domain of 
attraction of an SaS law. It was shown in IIKS79IIWan03IICD09ll that the following weak convergence 
holds inC([0,oo)): 

4? Znt inc. ,Wh')^ Af (Mo ,Y). (24) 

Henceforth we will use H' for the Hurst parameter of the collecting process and H for the Hurst 
parameter of the resulting stable process in random scenery. 

The Hurst exponent H = l — H' + H'/a can be explained by using the local time scaling relation 

(%(ci, x), a; e M, t > 0) = (^c^^^'iyit, x/c"'),x G M, t > . (25) 
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In IICS06L weak convergence in C([0, oo)) was shown for a properly normalized random reward 
scheme 

i=l 

where Zt is the linear interpolation of Z„ in the same manner as N are indepen- 

dent copies of mean zero, finite variance (/3 = 2) random walks which have H' — 1/2 explaining 
the exponent H = ^tl. They collect independent copies of i.i.d. rewards {ria\k)}i^jq which are 
also independent from the random walks. Cohen and Samorodnitsky called the limiting process an 
fBm-1/2 local time fractional stable motion. In tCD091, the discrete collecting process was gener- 
alized to Gh' and convergence to fBm-iJ' local time fractional stable motions for any H' E (0, 1) 
was proved. In IIDGP08L a collecting process scaling to /3-stable Levy motion (/3 > 1) was used, 
and consequently, other local time fractional stable motions were obtained in the limit. Let us now 
explicitly state these collective results. 

Recall that (fi', , P') is the probability space of Yf Suppose Mi{duj' , dx) is an SaS random 
measure that has control measure P' xLebesgue, but lives on some other probability space (fi, J^, P). 
As above, Wh' is either Gh' or Sp with /? e (1,2]. Letting H be as in ( |22] i, in light of ( pO] l we define 
a local time fractional stable motion as the process 

rf{Mi,Y)-^l eY{t,x;uj')Mi{duj',dx), t>0. (26) 
Jn'xR 

Let (c„) be an integer sequence with c„ — >^ oo and let {?7q ' (fc)} be independent copies of i.i.d. re- 
wards in the domain of attraction of an SaS law. The following weak convergence holds in C([0, oo)) 
as n — cx) : 

E ^ZrM^\W^^}) ^ rf (Ml, y) (indpt. seen.). (27) 

1=1 

Let M2 be a stable random measure with Lebesgue control measure with the restriction that 
a e (1, 2], and again let H be as in ( |22] i. We may use pO] ) and Holder's inequality to define 

Af(M2,f):= [ -E'lY{t,x;io')M2{dx), t>0. (28) 

JR 

Note that the scale parameter at time t for ( [28] ) is 

a= ||E'f^(t,a:;c^')llL°(K) (29) 

versus a ~ x; i^')\\L''{n'xR) for pS] ). For a £ (1, 2], a convergence result (in finite dimen- 

sional distributions) with respect to the single scenery case was given in Theorem 4.2 of IIDGP09I : 

E ^Z,,t{fji'^,wi^}) ^ Af (M2, y) (single seen.). (30) 

1=1 

As stated earlier, the process on the right-side is iJ-ss, but using (|29| one can see that this process 
does not in general have stationary increments. 

It is convenient to write p6l ) and ( |28| ) as renormalized sums of ( |23| l which appeal to the Stable 
Central Limit Theorem and the Law of Large Numbers, respectively (see IICD09I iDGPOSl IDGP09II ). 
The former renormahzation is applied to the entire integral in ( |23| ) and the convergence is inC([0, 00)) 
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whereas the latter renormahzation applies only to the integral kernel: 

n 

n-V"^(Af)« ^ rf, (31) 

1=1 

ll^n-'p^£f{t,x)^ Ahidx) m Af. (32) 

2.3 Scaling limits of alternating random reward schema 

We are now ready to state our results concerning the scaling limits of A{t) and its associated random 
reward schema ( [T7] i and ([T8|. 

Throughout this subsection we assume that the discrete collecting process Wh' [n) is extended 
to continuous time by linear interpolation and that it has the scaling limit Yt as given in condition 
(SLP). Independent copies of i.i.d. rewards {ria\k)}i^fi are, as usual, in the domain of attraction of 
an SaS law (scale parameter a — 1) and independent from the random walks. The space [Vl' . F' , P') 
supports Yt, and the SaS random measures Mi are as in the previous subsection. Define the processes 



AhW = AH(i;Mo,i^) := / l[o,y,(„o](2^)^o(rfa:), i >0, (33) 
TH{t)=TH{t;Ah,Y) := [ l[o,y,(^,)](x)A/i(dw', rfx), t > 0, (34) 
AH{t)^AHit;M2,Y) := ( Y.'l^o,Y,(^')]{x) M2{dx), t (35) 



which are analogous to (23 1, p6| ), and ( |28| l. 

The above are all self-similar with common index H = H' /a, and ( (34| l and ( (35| l are SaS pro- 
cesses. One can also observe (see Theorem 2.2 in llJunl II ) that both ([33| and ( |34j l have stationary 
increments. We call ( |33| l an (H'-sssi)-time a-stable Levy motion or more generally a stable process 
at random time. If X{t) is a two-sided a-stable Levy motion, then we may also write ( (33] l as X{Yt). 
The process ( |34[ i is an indicator fractional stable motion as introduced in IJunl 1 J . The process ( |35| l 
is the alternating analog of the scaling limit of a single scenery random reward scheme introduced in 
IIDGP09I . 

Theorem 2.3. Let H ~ H' /a, and let c„ — ?> oo n — >■ oo. 

• The following convergence holds inf.d.d.: 

n-"S^{WH'{nt)) Aff(i;Mo,y). (36) 

If the rewards are symmetric with finite variance (a = 2) and WH'{nt) converges in 

2?([0, oo)) ('C([0, oo))), then ^36\ also holds in T>{[0, oo)) ('C([0, oo)) respectively). 

• Ifn~^\WH{[nt\)\ is uniformly integrable, then 

Cn'/"£""''^"t(?7i^\W^i;,^) ^ T Hit-, M^,Y) (indpt. seen.). (37) 

i=l 

• Ifa>\, then 

Cn'£""''^nt('7a^W^i^'') ^ Ani^M^^Y) (single seen.). (38) 



3 PROOFS OF WEAK CONVERGENCE 



11 



The interest of the first convergence result (to (i/'-sssi)-time a-stable Levy motion) lies in the 
fact that this seems to be the first such Donsker-type theorem for iterated processes where the random 
time process is not a subordinator, i.e., not an increasing Levy process. In the case where the random 
time process is a subordinator, similar convergence results are well known. In fact, in Section 2.2 
of MNWXll I, such results are extended to the case where the scenery have a certain dependence 
structure. Their Donsker-type theorem shows convergence to an a-time fractional Brownian motion. 

It is not hard to see that A//(t), r//(i), and Anit) are all continuous in probability. However, 
by Theorem 10.3.1 in |ST94|, when a < 2, Anit) and F/f (t) are not sample continuous. In those 
cases, the best we can hope for is weak convergence in I?([0, oo)). We will see in the remark at 
the end of Subsection |3.2[ that even this is a lot to ask. In that remark, it is argued that even in the 
simplest cases A/f (i) is not even in I?([0, oo)). In particular, the weak convergence in C([0, oo)) and 
2?([0; oo)) given in the first part of Theorem |2 . 3 1 depends heavily on the fact that a = 2. In this case, 
the scaling limit of Sa is continuous since it is simply Brownian motion. 

The condition that n^^ Wh' ( \nt\ ) is uniformly integrable holds when Wh' is either Gh' or Sp, 
/3 > 1. The former follows from a Gaussian concentration inequality which bounds Wh' ( \nt\ ) 
in LP for all p > 1 (see fCM], p. 60), and the latter follows from eq. (5.s) in IIGR91I and the bound 
E(|X|U) < |lX||p(P(A))V«. 

3 Proofs of weak convergence 
3.1 Finite dimensional distributions 

A convenient tool in proving convergence of the finite dimensional distributions is a diagonal con- 
vergence theorem of |Doml 1 1. In order to state this theorem we require some definitions. 

As usual rja{k) is in the domain of attraction of the SaS law with scale parameter a — 1, and 
it is the reward on the edge between k and fc + 1. For fixed positive h, define to be the random 
signed measure on M which is a.s. absolutely continuous with respect to Lebesgue measure and 
whose random density is given by 



For < a < 1, we will say that (/„)„gN converges to / in 2?" if the following two conditions hold: 

• for any compact K C M, converges to fl^ in L^(M), 

• there is some ?/ > such that fn{x) = o(|a;|~'') and f{x) — o(|a;|^'') as a; oo. 

Let — L"(M) if 1 < a < 2 and J^" — 2?" if < a < 1. The following diagonal convergence 
is shown in Proposition 3.1 of IDoml II (see also Proposition 3.1 of [IDGP09I ). 

Proposition 3.1 (Dombry). Suppose Mo{dx) is an a-stable random measure, a € (0,2], and 
{fn)n<£N converges to f in T°^. If hn — > fli n — >■ oo then the random variables con- 
verge weakly as n oo and in particular, 




(39) 



For a locally integrable function / e L\^^, define 




(40) 




(41) 
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We now start by showing convergence in f.d.d. for Theorem 2.3 However, to reduce notation 
and simpHfy the presentation, we only prove convergence of the one-dimensional distributions for 
some fixed t > 0. The extension to f.d.d. in all three cases follows easily using the Cramer- Wold 
device (see for example Theorem 3.9.5 in HDurlOl ). 

Also without loss of generality we use A^^t] instead of the linear interpolation 
since they differ by at most n^^r]a{k) which goes a.s. to as n — > oo. 



Convergence in f.d.d. for ([36]l: Fix i e [0, oo). Let X„(i) = ■:^WH'{\nt\). According 
to assumption (SLP), X„(i) y{t)- By Skorohod's representation theorem, there is a common 

probability space on which X„ = Xnit), Y = Y{t) live and such that Xn{oj) Y{uj) for all 
tj e (note that the bar includes the dependence on t). 

Fix an ui and recall that H = H' /a and that for a < we let [0, a] :— [a, 0]. We have 



fcez 



(fc+i)i^ 



l[o,x„(<i)](2^) dx 



( [ntj )>fc>0}U{ W„/ ( [ntj )<fc<0} 



= n ^AL„tj(?7a,M^ff'(w)). 



(42) 



By Proposition 3.1 and the fact that Ijq x„{ui)] ^ ^[o y(o)] ™ ii3.ve that the one-dimensional 

distributions of A^nt] iVa, Wh') converge to those of Anit; A/q, Y). 



Convergence in f.d.d. for (38 1: For multiple independent walkers in the same scenery, we 
follow the arguments of Proposition 2.4 in IIDGP09I . Fix t e [0, oo). As in the proof of ( |36l ), using 
Skorohod's representation theorem and Proposition 3.1 we have for a £ (1, 2]: 



where Xn"^ (lu) — > Y'^'^ (cj) for each i G N and for alluj € Q (the bar includes the dependence on t). 
We need only show the following converges in probability to zero as n ^ oo: 



^ (l[0,x(-'(u;)]) ~ l[0,y(')(i;)] + l[O.Y(»)(a;)] - E'l 



(43) 



where for fixed n, the random variables Xn'' ,1 < i < Cn are i.i.d. Also, for each fixed i, Xn^ 
converges a.s. to Y^'^\ We first show that as n — > oo 

i—l 
1 II 



(")] 



(44) 



Consider a triangular array such that for each fixed n, there are c„ i.i.d. random variables 



<i<c„ 
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in each row, and for each fixed i, the column of random variables (C/j^"^)„gN converges weakly to 
zero. For such triangular arrays, the following weak law holds (see Proposition 2.4 in [,DGP09J ): 



as ri — > CXI, 



(45) 



thus proving ( |44] l. 

Since E||l[o.yj] ||La(R) < oo, the Strong Law of Large Numbers for Banach space valued random 
variables implies that the following converges a.s. in 



(46) 



thus completing the proof of one-dimensional weak convergence for ([38]). 



Convergence in f.d.d. for (|37]) : We will mimic the arguments of IIKS79I IDGP08I ICD09II . Let 

l{lV„,(nt;fc)}(w') := '^-{Wh' {[nti)>k>a}U{WH> {[nti)<k<0}i^') ■ 

Using the last equality in (|42]l, we have 



= Eexp £z0c;:i/"n-^^ryW(fc)l^^(., 



E' 



fcez 



}{nt-k)} 



where is the real-valued characteristic function of a symmetric reward rja and 



(47) 



(48) 



Suppose (paiv) ~ exp(— is the characteristic function of the SaS law of scale parameter 
(7 = 1. We show that the following asymptotic holds as n — oo: 



E' 



Lkez 



E' 



Y[ <paiVn{uj\ k)) 



Lk£l 



+ o(c„i). 



(49) 



If {xi)i^z and (x^iez are sequences in [—1, 1] with only finitely many terms not equal to one, then 



n-: n 



Letting 



g{y)^ sup \x\ "\(j),^^{x) - (pa{x)\ , x ^ 0, 

\x\<y 



(50) 



(51) 
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we have 



H ^Vc. {Vn{uj', k)) ~ 0a (Wn(w', k)) 

kei fcez 

< C„ ^ I k)) - Cj)a {Vn{uj', k)) \ 

feez 

< g{snp\vniuj',k)\)^Cn\vn{uj',k)\°' 



)^\n "eii^w„,{7U;k)}{^')\ 



kei 



(52) 



= gi9c-'^'^n-")\0rn-"'\WH'ilnt\;u')\. 

By assumption, rt^-^ | VF//' ([nij ;oj')\ converges weakly and is bounded in L^, so to prove ( |49| l 
we need only show that g{9cn^^"n^^) is bounded and converges in probability to 0. Since rja is in 
the domain of attraction of the SaS law with cr = 1, by the Stable Central Limit Theorem we have 
that as w — > 

<f>r,Av)^Mv)+o{\vn. 

Thus g is bounded, continuous, and vanishes at 0. Eq. (|49l) follows since 9cn ^^°'n^^ goes to zero. 



Let {{Ca\k))kez}eeN be independent copies of i.i.d. rewards such that Ca\^) has an SaS law 
with scale parameter a — 1. Using (|49|, the c„th root of ( |47j i is equal to 



Eexp (*0c-i/"n-^^Ln*J M^ff' )) 
= Ecxp (z0c-i/"„-^Al„*j (e^, <V)) + o(c-i) 

= E'cxp(-c-ln-«'^(^^l{H^^,(„,^fc)}(a.'))") 
= E'cxp(-c-in-^'|Ty^,(LntJ;c.')|0") + o{c-') 



(53) 



If bn is such that c„5„ — > A then (1 + 6„)'^ 



Letting bn = -c-iE'(n-^'|W^ff(LniJ))|0" 



and using the assumption of uniform integrability we have that 



(54) 



as required. 



3.2 Tightness 

In this section we will prove tightness for the family of processes on the left side of ( |36l ), thereby 
completing the proof of Theorem |2.3| Let us note that tightness results for iterated processes where 
the random time process is nondecreasing and converges to a stable subordinator are quite standard 
(see Chapter 13 in [WhiQ2|). The difficulty we face in this subsection is caused by the fact that our 
subordination is non-monotonic. 
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Proof of tightness for ([36]) when q = 2: We prove the case when Wh' (nt) converges in 
P([0, oo)); the C[0, oo)) follows similarly. Since we are in the case a = 2 we drop the subscript a 
and write 3 = 8^. Recalling H = H' /2, set 

X„{t)^n-''S{WH'{nt)). (55) 

By Theorem 15.2 in IIBil68l . it is necessary and sufficient to show 

(i) For each e > there is A/ > such that 

P( sup \Xn{t)\ > M) <e for all n and 

0<t<l 

(ii) For each e > there is an uq and < S < 1 such that 

P{m'x^ (S) > e) < e for all n > no- 
where m'^^ is the cadlag modulus. We note that in the C[0, oo)) case. Theorem 8.2 of MBil68l simply 
replaces m'^ in the above with the the continuity modulus mx„ ■ 

On the event {WH'{nt) e [m, m + 1)} we see that |X,i(i)| is bounded by n^"{\S{m)\ + \i]2\), 
the sum of m i.i.d. symmetric rewards plus an [772 1, a single reward which is independent of the sum. 
Choosing Mi such that 

P(h2| >Mi) <e/3, 
we may apply Levy's maximal inequality to get, for every M2 > 0, 



P ( sup \Xn{t)\ > A/i + M2 \WH'{n)\ £ [m, m + 1) ) < 2P{\S{m)\ > A'hn") + e/3. (56) 

\0<t<l 



Now choose C so that for all n 

00 

P(|W^ff'(n)| G Km+l)) <e/3, (57) 

which is possible since Wh' (n) satisfies (SLP). We then have, for each n 

P( sup \Xn{t)\ > M1+M2) 

0<t<l 

00 

< €/3 + 2Y,Pi\WH'{n)\ e [m,m + l))P{\S{m)\ > Man^) 

< 2e/3 + 2 J2 Pi\WH'in)\ e[m,m + l))P{\S{m)\> M2n"). (58) 

Next, since S{m) scales to Brownian motion, we may choose M2 such that P{\S{[Cn^ J)l > 
M2n^'/2) < e/6 for all n. We then have that (|58]l is bounded by 

2 ^'^""'^ 

j + l J2 P{\WH'{n)\e[m,m + l))<e. (59) 
Letting M = Mi + M2, we have proved (i). 
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We now show that (ii) is satisfied by Xn{t) = (nt)). Let n'~"'WH'{nt) 

and Sn{t) := S{nt)/^n. By convergence in I?([0, 1]) and (i), we may choose T so that 

P( sup |y„(t)| > T) < e/3 for all n. (60) 

te[o,i] 

By convergence in C([0, 1]) (see (ii)), given e > 0, one can choose 6i and rti e N such that 

P(m5^((5i) > e) < e/3 foralln> ni, (61) 

where the continuity modulus is taken over t e [— T, T]. 
Next, choose 82 and n2 G N such that 

P(m'y^ (Ja) > < e/3 for all n > 712, (62) 

where the cadlag modulus is taken over t E [0, 1]. 

Let {ti} be increasing times in [0, 1] satisfying ti — > 62- Setting m = , WH'{ns) — 
mu, and Wh' [nt) — mv, we have by (|60| and ( |62] i that for m large enough 



P {m'x {S2) > e) = P inf sup - S{WH'{ns))\ > e (63) 

\^{*Js,te[t._i,ti) J 

< 2e/3 + P sup ^— '-^^ ^ ^ 

\u,vel-T,T]:,\u-v\<5i V™ J 

By ( |6T| , the right side is bounded by e for m large enough. This proves tightness when a = 2. 

Remark: If a < 2, then Sa scales to an a-stable Levy motion, X{t). Fix e > and let > 
be the first positive time such that \X{Te) — \iui^_^^- X{t)\ > e. Consider the simple case where 
Wh' scales to a Brownian motion, Bf. Let t be the first time Bf — hits 0. As is well-known, 
Bt — Te oscillates around immediately, thus limj_j.^+ X{Bt) does not exist a.s. This argument, 
which can be made rigorous, shows that even in the elementary case where the collecting process 
scales to Brownian motion, the process X{Bt) is not cadlag. 



4 A recursive construction of some fractional stable motions 

Throughout this section we will suppose that a £ (1,2]. We present two related recursive construc- 
tions of some iJ-sssi processes. The first recursion produces stable processes in random scenery, 
while the second recursion produces local time and indicator fractional stable motions. Note that 
only the second recursion leads to SaS processes. Since fBm is the only sssi Gaussian process, when 
a = 2 the second construction gives us fBm. In particular, if on the first step of the recursion we 
use Brownian motion as the collecting process (or random time process), then we obtain fBm of any 
dyadic Hurst parameter 

Although the first construction does not in general lead to a-stable processes, we will see that 
the finite dimensional distributions of the processes have finite a moments, and thus one can appeal 
to the Stable Central Limit Theorem and normalize partial sums of independent copies of the stable 
processes in random scenery in order to get honest stable processes (in a mann er si milar to ( (3T] i) . 



Let be an i7-sssi process satisfying the four conditions of Theorem 



4.1 



below. Consider 



the vector v = {vi, . . . , Vn) with coordinates vj G {+, — }. Let us use the notation v to denote v 
truncated by removing the last element, that is, w = {vi, . . . , Vn-i)- The empty set will denote the 
empty vector 
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We define the process F^" recursively from Yf" and an a-stable random measure Mo{dx), with 
a e (1, 2], assumed to be independent from Y^^ . If w„ = (+) we let 

F," := f eY4t,x)Mo{dx), (64) 

and if Vn — {—) we let 

y,'' := / l^,^y.^{x)M„{dx). (65) 

The second recursive procedure is defined similarly. We again use vectors, now denoted w — 
{wi, . . . , Wn), with coordinates taking one of two different values. However, in order to distinguish 
between the two procedures, we let Wj G {*, x}. As before, we let w = {wi, . . . , Wn-i). 

Once again Y^"-' is defined recursively from Y^ and an a-stable random measure Mi with a G 
(1,2], however, the control measure of Mi is no longer Lebesgue measure as it was in the case of 
Mq. Suppose that {fl' , P') is the probability space of Y^*". Then, just as in (|26|, Mi {du' x dx) 
has control measure P'xLebesgue and lives on some other probability space (fi, P). If w„ = (*) 
we let 

yj"" := / iY^{t,x){uj')Mi{duj' X dx), (66) 

and if w„ = ( x ) we let 

Yr--=[ l[o.Yfi.')]i^)Mi{du;' xdx). (67) 

JSl'xR 

We must show that the above recursions makes sense, i.e. that the integrals are well-defined. In 
general, it is known that _ff-sssi SaS processes have local times almost surely. This almost 

gets us to where we want to be, however there are two separate issues with which we must deal. 

According to (19\ we need that the integral kernels of (|64]) and ( |65| l are in L" (M) (which easily 
follows if they are in i^(M)), but (|64]) and ( |65| l are not in general SaS processes, and thus we need 
an extra argument to show that they have L^(M) local times almost surely. 

The second issue concerns (|66| and ( |67] i which are SaS processes, but are well-defined only if 
the local times are in L"(ri' x M). In other words, we will need the a-th moment of the local times 
to be integrable. To solve these two issues, we use the following result. 

Theorem 4.1. Suppose Yt — Yt{Lo') is an H'-sssi process which satisfies 

(a) < E|Yi| < oo. 

(b) Yt has a local time satisfying < E J^f-Y [t, x)^dx < oo. 

(c) Yi has a bounded continuous density. 

(d) E[suptg[o^i] \Yt\] < oo. 
Then the processes 

= / eY{t,x)Mo{dx), (68) 
y}-^ = f ^o.Y]ix)Mo{dx), (69) 

^R 

y}*'^ ^ I £Y{t,x){uj')Mi{duj' X dx), (70) 
Jn'xR 

Yl''^=f l[o,Y(u')]{x)Mi{dL,' X dx), (71) 
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are well defined H-sssi processes satisfying (b) - (d), where H = 1 — H' + H'/afor and 
y(*^ and H — H'/afor y(^) and Y''^\ Moreover, all four processes have finite a moments which 
implies they also satisfy (a). 

Remarks: 

1 . For the proof, we need that f y satisfies the occupation time formula 

lA{Ys)ds= I eY{t,x)dx, 

J A 

for any Borel set A. This follows from definition ( [2T] i. 

2. The processes y/^'' and 1"/ ' are not in general stable. However, as mentioned above, when 
they have finite a moments, one can use the Stable Central Limit Theorem and normalize 
partial sums of independent copies of these processes to get stable processes. 

Proof: 

Well-defined i7-sssi processes with finite ath moments: To see that the Y* are well defined 
and satisfy (a), we have 

E[{YtY] = E[(f/" V] = E' / \l[o.n]{xrdx = E'lni, (72) 

Jr 

which is positive and finite since Yt satisfies (a). Also, 

E[(r/+V] ^ E[(y/*V] =-E' J £y(l,a;)"dx. (73) 

To see that ( |73] l is finite and nonzero, note that E' (1, x)dx = 1 by the occupation time formula 
and!]' J^iy {I, xfdx < oo by (b), thus ly € i"(rj' x M) f or a e (1,2]. 

To see that the Y* are H-sssi, we refer the reader to Theorem 3.1 in IICS06II and Theorem 2.2 in 
OunTTl . 

Property (b): Next we use Theorem 21.9 of IIGH80II which implies condition (b) under the as- 
sumption that 

f f f E[e'^'^^' ^^'^]ds dt d9 <oo. (74) 
Jr Jo Jo 

(*) 

Let us show (b) for . We have 

E[e*«(^**"-^i*')] =exp (^-rE'y |^y(t, x) - (s, a;)rda;^ . (75) 

Using iy^t, X — Bs) — £y{s, x — Bs) ^ ly{t — s, x) and 



we see that 



ey{ct, c"'x) = c}-"'ly{t, x) (76) 



\ty{t,x)~ly{s,x)\°'dx = liy{t~S,xYdx 



\t - s|"(i-f^'+^'/")£^(i^ u^du. {11) 
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Substituting v = e-\t~ s|i--ff'+-ff7" (e' J^£y{1, we get that ([74} equals 



which is finite since E' (1, u)°'du > by the occupation time formula. 
To show (b) for y}^\ write 



(78) 



cxp (^ie J (eyit, x) - iris, x))Mo{dx)^ dujduj' 

E'exp(-r / \tY{t,x)-lY{s,x)\°'dx 
\ Jr J 



(79) 



Using ( |77] i and ( |79| l we have that in this case ( |74| i is 



E' 



exp -r / \t - s\°'^^- "'+"'/ °'kY{l,u)°'du] d0 



ds dt. 



Substituting v ~ 9 ■ \t — ^ u)°'duj " and integrating we obtain, for some 

constant c > 0, 



cE' 



-l/a 



T /.T 



dsdt 



Jo 



To show that this is finite we need only show that E' (Jjj i'y (1, 
^y(l,a;)(a;') = Ofor 

\x\>A{uj') := sup |rt(c^')L 



te[o,i] 



so by Holder's inequality 



A(u,') 



< 



£Y{l,x){uj')dx ^ / ^y(l,x)(a;')fia; 



(80) 

< CXI. We have 
(81) 

(82) 



By the occupation time formula, the left side of ( (82] l equals 1 a.s. so that 



E' (^J iYil,x)°'dx^ ' <E'(2yl)(' 



Q— 1) /a 



(83) 



Property (d) of Yi completes the proof of (b) for F^^^. 



Moving on to , we have 

= exp (-rE'|rt_,|) = exp (-^ |i - sfE!\Y^ 
Thus ( |74j l reduces to 

rj-i rj~i rj-} rp 

[ [ [ exp(~0°'\t-s\"'E'\Yi\)dtdsde = C [ [ \t- s\-"'/"dtds 
Jr Jo Jo ^ ' Jo Jo 



(84) 



< oo, (85) 
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where C = exp(-w"E'|Yi|) du. 

Finally, let us consider y (~) . We may mimic steps ( 75 1 through ( 80 1 in order to reduce ( 74 1 to 
showing 

-l/a . 



E' 



(86) 



But this follows from assumption (c) on Yt, since we may simply integrate |a;| against the 
bounded continuous density of Yi which will give a finite value. This establishes (b) for Y'-^\ 

Property (c): In the course of showing property (b) for Y*, we showed that in all cases Y* 
possesses a nonnegative and integrable characteristic function, and thus (c) follows from Theorem 
3.3.5 in IDurlOJ. 

Property (d): Consider first a = 2. Property (d) is known for y}*^ and Y^^^ since they are sssi 
Gaussian processes, i.e. fractional Brownian motions. 

For y}^\ let Bt be a two-sided Brownian motion. We use Proposition 2.2 in I1KL98II which is 
essentially a corollary of Slepian's Lemma. It implies that for each fixed uj' 



P sup / l[o^YAu:')]is)dBs > y ] < 2P / l[o,Y,(ui')]is) dB, > y 
\te[o,i]JR J \Js. 

Integrating over Q', property (d) for y}^^ follows from property (d) for Yt. 



(87) 



For Yt , let Y* := supti^in.i] ^t, and := inf(g[o.i] Yt- We have 



E[ sup \Yt 

tG[0,l] 



(-)l 



< E[ sup y/ ^ + sup (-r/ ')] 

te[o,i] te[oa] 



< 2E'[ sup / l[o,y,](s)dB,] 
te[o,i] 



< 2E'[ sup / lio,T]{s)dB,] 

Te[Y,,Y*' " 

< 8E'(r*). 



(88) 



The last inequality follows since the integral in the second to last line is just a two-sided Brownian 
motion at time T and E'(y*) = E'(— F*) < co. We thus get property (d) for 5^*'"^ since property 
(d) holds for Ft. 

Let us now suppose that 1 < a < 2. Theorem 10.5.1 of IIST94I states that if 



Yt - / ft{x)M{dx) 



(89) 



for some family of L°'{E, m) functions {ft{x)}t>o, where m is the control measure of M, then there 
is a constant C such that 



P( sup \Yt\ >y)<^ [ sup \MxTm{dx), 
te[o,i] y JEte[o,i] 



(90) 



for any y > 0. 

We can therefore obtain (d) for Y' by showing that 



E' [ sup (^y(t,x)")dx = E' [ lY{l,x)"'dx 



(91) 
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and 

E 



' / sup (l[o,y,](x)")dx = 2E' ( sup Yt] (92) 

JR«G[0,1] \t6[0,l] / 



are both finite. As seen in d73]) and the argument thereafter, (91 1 is finite since Yt satisfies (b). Also 



( [92| i is finite since Yt satisfies (d). □ 

For fixed a G (1,2], define 

(j>+{x) := 1 — X + x/a and (t)^{x) := x/a. (93) 
Applying Theorem |4. 1 [ recursively, we have the following corollary: 



Corollary 4.2. // is an H'-sssi process satisfying (aj-(d) of Theorem \4.l\ then y^^^i Qjid 
y(wi,...,w,i) defined in \6A) , \65) , are H-sssi processes with 

Moreover, - q„ SaS process. 

5 Brownian motion extracted from fBm, H < 1/2 

Suppose a — 2. Then the family of stochastic integrals, {Y}^^)t>Q, is an i/'-sssi Gaussian process, 
thus it is precisely fBm with Hurst exponent H' < 1/2. In this section, we show that Brownian 
motion can be extracted from Y^ by time-changing its integral kernels. In order to motivate our 
time-changed kernels, we first show that Brownian motion can also be extracted from a stable process 
at random time, y} \ using a time-change. 

To keep things simple, we assume in this section that the random time process Yt is itself an fBm. 
Thus it is a.s. continuous and satisfies the property that for each s > 0, 

Ts = infjt : Yf = sj < oo a.s. (94) 

Heuristically, time-changing the kernel of Y* undoes the subordination of Y* to the process Yt, 
leaving us with a process (M {At))t>o- We then observe that As C At for s < t, and that m{At) is 
linearly increasing (here m is the control measure). One need only check that such a procedure gives 
us what we want, by looking at the finite dimensional distributions. Since our interest is in the case 
a = 2, we have that Mq, Mi are Gaussian random measures on M and ft' x M, respectively, and we 
in fact need only check covariances. 

Let us start by presenting the time-change of Y^ ^ . 

Proposition 5.1. Let the random time process Yt be a fractional Brownian motion. IfY^ ^ is defined 
as in ( |69| l with a — 2, then Yt^ ^ is a Brownian motion. 

Proof. We have 

^ I l[a,Y^,]{x)Ma{dx) - / l[o,t](s)rfi3. = Bt, (95) 
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where Bt is a two-sided Brownian motion. For the covariances, if s < twe have 

(l[o,,](r))2dr = s. (96) 



□ 

In the case of 

Jn'xR 

we cannot look at "ly^ ^" since Tj hves on the same probability space as Mi. We address this issue 
by instead time-changing the kernel l[o,Ft]- Let us define 

l[o,y^^(„,)](x)A/i(d^' X dx). (97) 
A good way to think about the above integral is in terms of a Central Limit Theorem similar to ( (3T| : 

n „ 

n-i/2^Aff(T«)« ^ / l[o,^^^(^,)](a;)Mi(dw'xdx). (98) 
._i Jn'xm. 



Here, r^*) is measurable with respect to the cr-field of A^-*. By Proposition 5.1 the A/^ (r^ ) are 
independent Brownian motions. The next proposition shows that the right side is also a Brownian 
motion thus proving ( |98] l. 

fx) 

Proposition 5.2. Let the random time process Yt be a fractional Brownian motion. IfY^ is defined 
as in with a — 2, then l"/^"*^ is a Brownian motion. 

Proof. We have 

r/^)- - / l[o,Y^j{x) Ahidu' X dx) 
Jn'xR 

llo.t]{x) Miidco' X dx) = Mi{n' X [0,t]), (99) 

nxR 

which is a Gaussian random variable with variance P'xLeb.(i7' x [0,t]) — t. For the covariances 
we analyze second moments. If s < t we have 



(ri^)- + Yy>-y - E [ / (2 • l[o^,j + l[s,t]) MiidLo' X dx)] 

= / (2 • l[o,s] + Ms,t]f P' X Leb.(da;' x dx) 

Jn'xR 



E 



(4 • l[o,s] + l[s,t]) dx 
= 3s + t (100) 

as required. □ 
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